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ABSTRACT 


The development of new composite materials, which lack the historical field data base, has led 
to the need for an accelerated life testing method applicable to composites. Accelerated life testing 
by increasing the sustained stress level requires the modeling and validation of a strength-life relation. 
Proof testing of composite fibers by over-loading is one step in the understanding of the relationship. 
It is also important in the reliability and safety assurance in deployment of composite structures. 

A parametric study examined the strength life relation of composite fibers and a methodology 
to analyze the fiber failures after proof testing. The fiber statistical strength was modeled by a 
probability of failure model, while a deterministic approach was taken when considering individual 
fibers and the associated life reduction each fiber experienced during the proof testing procedure. 
Also studied was the distribution of the first failure to occur after proof testing in order to understand 


the effects of the sustained load and the proof load on fiber life. 
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I. INTRODUCTION 

Aircraft structures built of composites materials are 
becoming the norm rather than the exception today. As 
stronger fibers are developed in conjunction with new matrix 
materials, the advantage of higher strength to weight and 
stiffness to weight ratios become more apparent. The dynamic 
development of composites necessitates a predictive 
methodology for safety and reliability in the absence of 
historical data. A designer using a homogeneous material such 
as steel or aluminum only has the option to vary geometry in 
order to meet the required loading on the structure. Using a 
composite material, not only can the geometry be varied, but 
the material itself becomes a part of the design process. The 
design of composites must include a parametric understanding 
of the strength to life durability relation. Unfortunately, 
the historical data regarding safety and reliability is not 
available for newly developed composites. 

Mathematical modeling of the strength life relation must 
be based on the failure process of the specific material. The 
failure process for composites is inhomogeneous, typified by 
local fiber failures at isolated sites. The load sharing 
ability of the matrix binder will temporarily delay 
catastrophic failure of the composite. As the stress or time 


increases, the number of failure sites increases and chance 





clustering of such sites ultimately leads to failure. 
Experience in composite strength modeling has demonstrated 
that given the fiber strength statistics, the probability of 
failure of the composite can be predicted. Since the fiber 
failure process is homogeneous, the modeling of inhomogeneous 
composite life will be based on the fiber statistics. 

Currently being conducted at the Naval Postgraduate 
School’s Advanced Composites Laboratory is a strength-life 
experiment involving two sets of statistically identical AS-4 
graphite fibers. The life statistics of the fibers, which 
have been subjected to a sustained tensile load, are being 
collected and analyzed. Concurrent with the life testing, a 
proof test was conducted on two subsets of the fibers in order 
to validate the strength life model and to allow for 
accelerated life testing. 

The focus of this study is the life sensitivity to 
parameter variation and the post-proof test minimum time to 


failure life distribution. 





II. BACKGROUND 


A. APPLICATION TO MISSION 

1. Structural Safety 
Structural safety may be considered as the number of 
failures which occur during a given time period. For example, 
the failure of a composites helicopter blade may be modeled by 
an appropriate probability density function (pdf) as shown in 


Figure 1. 


Probability Density Function 
(pdf) 
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Pigure 1. Structural Safety 








The risk associated with a specific blade is then the 
probability that a blade which has lasted a certain number of 
flights, may fail during its subsequent flight. Therefore the 
height of the curve is a measure of safety. A low height 
means a safer blade and conversely, a high region indicates a 
relatively unsafe blade. 

The time period used to define the safety of a blade 
is one flight or a block time, but no two flights produce the 
same stress history on a blade. Therefore the block time 
must be convoluted to a reduced time at an equivalent stress. 
The idea of a reduced time was introduced by Coleman [Ref. 1] 


as 
Ce 
jel 
= (°) z{RiSKk)) at 


That is, given a time dependent stress history, S({), and the 
physical breakdown process, K, the fractional life of the 
fiber used between t, and ty is 7, where € is an intrinsic 
normalizing time parameter. An example of reduced time fora 
stress history is given in Appendix A. 
2. Availability 
While the probability density function (pdf) is useful 


in estimating safety, it can only provide guidance to 
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Pigure 2. Availability 
logistics and planning associated with the availability of a 
fleet of aircraft. If the pdf for the set of helicopter 
blades is known, then the cumulative density function (CDF) 
may be found by integrating the pdf (Ref. 2]. The CDF in this 
example is then the total number of blades, F(t), out of the 
original set that have failed up to a given time. Also, 1- 
F(t) is the number of blades available, R(t), at any time (see 
Figure 2). Early on there are many blades available but as 
the blades near the end of their planned service life, the 
inventory has dropped. The availability curve may then be used 
as an acquisition planning tool in order to optimize purchases 
of new equipment. 
3. Feasibility 

Another utility of the pdf can be found in estimating 

the success of a mission, or a measure of the feasibility. 


The instantaneous failure rate of the blades, A(t),may be 











dividing the pdf (# of failed) by R(t) (# available) as shown 
in Figure 3. The success rate, SR(t), can then be defined a 
1-A(t). The success rate may be used in conjunction with the 
availability curve in order to determine mission feasibility. 
For example, a certain mission that requires 100 aircraft 
would need 100/SR(t) aircraft to be deployed in order to 
ensure mission accomplishment. Looking at the number of 
aircraft available at that time would tell if the mission were 
possible. As the failure rate increased, there would come a 
point where the number of aircraft required would exceed the 
number available. By knowing the pdf, one could then predict 
not only the safety associated with a structure, but when a 


replacement must be brought on-line. 


Instananeous Failure Rate 
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Figure 3. Feasibility 








B. BASIS FOR PROOF TESTING 
1. Homogeneous Brittle Material 

Given that a material’s durability (life under a 
stress history) inherently contains statistical variability, 
there always exists a finite probability of failure. The 
probability of failure depends on the location of the time 
period with respect to the mean age. The underlying idea of 
proof test is to overload each sample so that the weak samples 
will be failed during the overload and therefore eliminated in 
the subsequent deployment. This assumes that the samples 
which survived the proof testing are not significantly damaged 
by the overload (in strength) and therefore its effect on the 
life of the sample is minimal. A rational proof test 
methodology must be able to quantitatively characterize the 
damage during proof test. Such methodology cannot be 
empirical because one sample cannot be tested in strength and 
then in life. An understanding of the failure mechanism is 
needed. 

For a homogeneous brittle material, the failure 
mechanism is flaw growth. Given a specimen with a crack of 
length a, the strength of that specimen can be related to the 


crack length by classical fracture mechanics as 


k=o/a 


where k is the stress intensity factor. A specimen with a 
short crack is stronger than a specimen with a longer crack 
for the same stress intensity factor (see Figure 4). At any 


time, the length of the crack may be determined by 


[- | sample 1 
[=] sample 2 


service stress 





Figure 4. Homogeneous Crack Growth 


where 6a/6 is the rate of crack growth. But the rate of 
crack growth is itself a function of the applied stress oa, 
i.e., 6a/6§ = alo). 

Considering two specimens as shown in Figure 5, it is 
desired to apply an intermediate proof load at ty, which would 
cause the weaker sample to fail but not the stronger. By 
failing the weaker specimen in a controlled environment, it is 


removed from service and prevented from failing while in 





service stress 





Figure 5. Proof Test of Homogeneous Material 


operational use. If the breakdown function is known, the 
absolute safety up to time t, may then be guaranteed. 

If the crack growth rate for each sample is the same 
and the stress level is identical, i.e., a(o), = a(o)., it 
implies a homologous relationship between the samples. Under 
the condition of homologous, the effect of proof testing may 
be examined in a deterministic setting. Physically this 
implies that each sample with an intrinsic (but different) 
strength will have a strength degradation contour as depicted 
in Figure 5. Given this hypothesis, a high proof load at an 
early time t, will assure safety up to time t,. This proof 
test strategy has two deficiencies. One, the high proof load 
will eliminate a larger number of samples, which may otherwise 
provide some useful service between t, and t,; a deficiency of 


economics. Two, the homologous condition may not hold for an 


extended time period; thus, the safety up to tz is not 
guaranteed; a deficiency in efficiency. 

The proof test methodology may be modified to include 
several strategically placed proof loads of smaller magnitudes 
dispersed through out the service life of the structure 
(Figure 6.) This method would involve a shorter time 
extrapolation, allowing for piece-wise homology and some 
useful service would be obtained by samples that otherwise 


would have been failed by one large proof load. 


: service stress 





Figure 6. Multiple Proof Tests 


2. Heterogeneous Material 
The failure mechanism for a composite structure is not 
homogeneous. Rosen ‘Ref. 3) has modeled the load transfer 
between fiber bundles and matrix by the Local Load Sharing 
model. The fibers carry the stress applied to the structure 


while the matrix transfers the load to adjacent fibers. When 


10 





a fiber fails, the longitudinal stress at the break drops to 
zero while the shear stress reaches a maximum as characterized 


in Figure 7. At the fiber break there is a region known as 





Figure 7. Local Load Sharing 


the ineffective iength, 6, over which the fiber does not carry 


any longitudinal stress. The ineffective length is given by 


2-4 )2 -1 
8=((V, uz cosh (2etare) >) ld, 


where: Vz, is the volume fraction of the fiber 
E, is the modulus of the fiber 
G, is the shear modulus of the matrix 
o- is the fractional value, called the fiber load 
sharing efficiency, below which the fiber is 
considered not effective. 


11 








If enough fibers fail in close proximity to each other, this 
clustering of failures can lead to catastrophic failure of the 
composite. 

The composite strength can be modeled by knowing the 
fiber strength, matrix ability to transfer load and the 
fiber/matrix interface. Harlow and Phoenix [Ref. 4 and 5] 
have developed a model which predicts the probability of 
failure of a composite in tension based on fiber strength 
statistics and ineffective length. The model has been 
extended by Harlow and Wu [Ref. 6] to include multi-modal 
fiber failure model, which has been verified by Wu [Ref. 7], 
Storch [Ref. 8], Kunkle [Ref. 9], Englebert [Ref. 10] and 
Johnson [Ref. 11]. 

In the case of a composite helicopter blade, a proof 
load applied directly will lead to unnecessary failure of 
fibers because of the load sharing. Any over load would lead 
to clustering of failure sites in the structure as the stress 
on surrounding fibers increased. By removing or reducing the 
load transferring ability of the matrix during proof test, the 
clustering would be avoided. A possible method to reduced the 
influence of the matrix would be by heating the composite 
structure (Wu [Ref. 12]). 

3. Composite Fiber 


B. Coleman’s failure potential theory is based on the 


*weakeat link" idea and is a stochastic function between fiber 








strength and fiber life (Ref. 1]. If an individual fiber is 
considered, the fiber has either failed or remains intact, and 
the failure is described by a binomial distribution. Fora 
fiber material with a high modulus, the failure mechanism is 
flaw growth. If the flaw is greater or equal to a,,;, the 
fiber fails, and if the flaw is less than a,,;,, the fiber will 
survive. As the diameter of the fiber becomes smaller, so 
does the allowable flaw size. If the fiber is divided into 
Many equal volumes, referred to as the metric, and the flaw 
density within each metric is low, the binomial distribution 
becomes Poisson. The probability of failure at t=t, , given 
a stress o, is F(t=t,|¢) = 1-exp(-¥(7)). Because the failure 
of the fiber is homogeneous to mechanism (flaw growth) the 
hazard ¥(r) is of the Weibull form. That is, ¥(7) = 7%, where 


Tt is the reduced time as stated earlier. 
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III. EXPERIMENTAL PROOF TEST PROGRAM 


Based on the background summarized, an experimental 
program is designed to collect graphite filament life data 
under sustained constant loads. The purpose of the 
experimental program is three fold. One is to characterize, by 
actual data, a formulation of the breakdown rule K. For 
example, if K is of the power law form, what are the values 
for t and p. Two, to examine the validity and range of the 
strength-life homologous correspondence. And three, to 
explore the effect on life of the proof test overload. 

A total of 512 test stations for filament life testing 
have been designed and constructed. These stations can apply 
a sustained load to a single filament fiber sample using dead 
weights with provision to isolate external disturbances such 
as those caused by earthquakes. Fiber filaments from two 
spools of graphite fibers (AS-4 manufactured by Hercules 
Corporation) have been put on sustained load, Ref. 12. These 
two spools were selected because they have different 
statistical characteristics in strength, as noted by Englebert 
[Ref.10]. Allocation of the test stations are equally divided 
into 256 stations for each of the two spools. Each of the 256 
stations are further subdivided into three load levels as 
indicated in Figure 8. Some of the fibers at the highest load 


have been realized in life while most of the samples at the 
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Figure 8. Fiber Life Testing Load Levels 


lower load levels are still surviving. In time, with an 
adequate number of failures as depicted by the data band, the 
parameters for the breakdown rule K can be assessed. 
Simultaneously, the two spools with different strength 
characteristics may produce two different data bands as 
depicted in Figure 9. Comparison of these two data bands by 
appropriate probability based statistical methods can confirm 
or reject the homologous damage hypothesis, i.e., longer life 
is associated with higher strength. 

Finally, in this investigation, proof test by overload is 
performed on half of the surviving samples under the highest 
sustained load level. The subsequent life of the samples 
surviving the proof load can then be compared to those samples 
that were not subjected to the proof test as illustrated in 


Figure 10. 
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Figure 10. Post-Proof Test Life Distribution 





Test stations vacated from proof testing are re-utilized 
by putting new samples on at the lowest load level to extend 


the prediction of the damage function. 
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The experimental procedure is described in Appendix B and 
the computer software used during the proof testing is listed 
in Appendix C. 

Up to now, sufficient time has not yet elapsed for the 
observation of samples failed from either history. In the 
future, as the data becomes available, they may be 
progressively interpreted using the analytical techniques 
described in the following parametric studies, replacing 


simulated data with actual experimental data. 
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IV. DETERMINISTIC PARAMETRIC LIFE STUDY 


A. FRACTIONAL LIFE 

Despite the large replication of samples being tested in 
the experimental program, relatively few failures have been 
realized in time; i.e., failures classified as life data. 
This is because of the interaction between the strength 
scatter and the insensitivity of the breakdown rule, K, on the 
sustained stress level as shown in Figure 11. The high 
strength scatter can be attributed to the small diameter 
(approximately 6 microns) of the fiber filaments in that even 
a sub-micron imperfection represents a large percentage of the 
load carrying fiber cross section. The insensitivity of the 
breakdown rule to stress level can be attributed to inherent 
stable carbon micro-structure; that is, the addition of time 
does not contribute to the kinetics of micro-damage. These 
two characteristics makes experimental planning difficult. If 
a high sustained stress level is selected, a large portion of 
the samples will fail during the loading process (realized in 
strength). If a lower sustained: stress level is chosen, 
because of the shallow slope of the strength degradation, an 
even longer life (logarithmic) is expected. The total effect 
is that mathematical and experimental ingenuity cannot reduce 


the time before the realization of data. 
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Figure 11. Strength Life Distribution Relation 


In this section, a numerical parametric study is made on 
simulated data. The underlying function used in the Monte 
Carlo simulation is a two parameter, Weibull, weakest link 
distribution for fiber strength based on the justification in 


section II.B.3, 


F(9;a,B) =1-exp(-(5)°) 


The numerical values of a = 5 and 8 = 20 grams are based on 
previous test data by Englebert [Ref.10]. 


The breakdown rule is assumed to be of the power form, 


~((t) 
Ra) 
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and the justification is the stability of carbon minimizes all 
kinetic damage processes (the exponential rule). The slope of 
the power law, p = 40 is selected from experience obtained 
from other fibers. Other parameters used in the simulation 
such as loading rates are those used in actual testing. 

Three stress histories were used, the first was a constant 
loading rate where the fibers were loaded until failure 
(Figure 12). Knowing the load rate and the failure strength, 
the failure time was calculated (Appendix D). 

The second stress history considered was a constant 
loading rate and then maintained at a constant stress level 
(Figure 13). By calculating the fractional life used up in 
each portion of the stress history, the failure time was 


determined (Appendix E). 
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Figure 12. Constant Loading Rate History 
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Figure 13. Stress Rupture History 


Finally, a proof test history (Figure 14) was looked at 
and failure times were computed. The notation used in Figure 
14 is defined in Appendix E. 

As a fiber passes through each region in Figure 14, a 
portion of its life is consumed. Although each fiber 
experienced an identical stress history up to failure, the 
fractional life used up in each region varies greatly because 
of the dependence the intrinsic strength of the fiber. A weak 


fiber uses a large portion of its life just to attain the 
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Pigure 14. Proof Test History 
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sustained streas level. But the initial loading is small in 
proportion to the intrinsic strength of a strong fiber, 
thereby having little effect. Figures 15 and 16 show the 
fractional life consumed for the same three _ fibers, 
representing a low, medium and high intrinsic strength, during 
two different proof test histories. As seen in both figures, 
the weak fiber’s entire life was consumed in just the initial 
loading process. While in contrast, the stronger fiber in the 
front of the graph, used a negligible amount of its intrinsic 
life reaching the proof load. By increasing the loading and 


unloading rate of the proof test portion from 0.8 and -0.8, in 


etre 08 


ody 





Figure 15. Fractional Life Run # 1 
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Pigure 16. Fractional Life Run #2 


run #1 to 50 and -50 (qm/sec), the time spent by the fiber in 
that region wan decreased and the fractional life consumed was 
also reduced thereby allowing the fiber to last longer in real 
time. Ry manipulating the parameters, the fractional life 
consumed in any portion could be controlled in such a manner 
as to selectively screen out the weaker fibers. Appendix F 
contains a listing of the software used to calculate the 
fractional life along witha listing of the parameters used in 


each run. 
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V. POST-PROOF TEST SAFETY ZONE 

As discussed in section III, for a period following the 
proof test, no failures occurred as compared to the non-proof 
test history. By causing the fibers to fail during proof 
testing, the time immediately following the proof test becomes 
a safety zone. A computer simulation was conducted to study 
the distribution of minimum times to first failure after the 
proof test, under varying parameters. A total of nine cases 
were conducted, each with 50 runs and each run simlated 2000 
fibers. A histogram was then produced for each case in order 
to determine the sensitivity to parameter variation. The 
software used to conduct the simulation is listed in Appendix 
G along with a listing of the parameters for each run. 

Figure 17 shows typical failure times for the non-proof 
tested and proof tested histories. Increasing the magnitude 
of the proof load for each sustained stress level caused the 
time to next failure to increase. The sustained stress level 
also played a role by causing the location of the histogram to 
shift. At lower sustained stress levels, the time to next 
failur= was larger than at the higher loads. Figure 18 shows 
the histograms developed using a low sustained stress level 
and three levels of proof load; high, medium and low. As 
expected, by increasing the proof load, the time to next 


failure increased. The effect of screening out the weak 
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fibers, those that would have failed in time, increases the 
safety zone tremendously. Figures 19 and 20 show similar 
trends for higher sustained stress levels. A high sustained 
stress level and a low proof load resulted in the minimum time 
before failure due to much of each fiber’s fractional life 
being consumed just to reach the time the proof test was 
conducted. This in conjunction with a non-rigorous proof 
test, which did not screen out the weaker fibers, yielded the 


shorter times to next failure. 
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Pigure 17. Stress History Failure Times 
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Figure 18. Post-Proof Test Safety Distribution Stress Level 
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Figure 19. fonrt-Proof Test Safety Distribution Stress Level 
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Figure 20. Post-Proof Test Safety Distribution Stress Level 
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VI. CONCLUSIONS AND RECOMMENDATIONS 
The process of proof testing a material in order to ensure 
that it will not fail when in use at its service stress level 
is a much more complicated procedure than merely choosing an 
appropriate proof load. By understanding the role of each 
parameter in the strength-life relationship, a proof test 
method can be developed which may yield very precise results. 
A specific recommendation is to include using a 
probabilistic approach with a model based on the pending data 


developed in this investigation. 








APPENDIX A. REDUCED TIME CALCULATION 
Given two stress histories as shown in Figure 21, what is 
the effect of sequencing of the stress levels and what effect 


does the slope have on the equivalent life? 


stress stress 





Figure 21. Stress History 
The reduced time given by 


74 
jel 
T (+) g [R(S(E)) a 


and is a measure of the intrinsic life of a specimen. K(®*) is 
the breakdown rule and S(é)) is the stress history that the 


specimen has undergone. 
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For this example, choose K(*) of the power form. 


y= (S08) 
K(+)=(=3°)e 


For case 1: 
S(g) = 12 for 0 s & < 50 
= 8 for 50 = € «< 100 
1p 12 ars 8 
-)z sf \P —)?P 
T(-) +f 2) doz fp) dé 
and evaluating at the limits of integration gives 
+) zt 
T(-) eye (12° (50) +8°(50)) 
For case 2: 
S(&) = 8 for 0 «= ~ < 50 
= 12 for 50 s € «< 100 
int 8 Phas 12 
aye &ye =e he 
T(+) 4p aoe fp) dé 
and evaluating case two at the limits of integrations gives 


1 
T(:)= Fa 8? (50) +12°(50)} 





Therefore, sequence is not a factor. 
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Now consider case 1 to find the reduced time at an equivalent 


stress S,. 


1 
T(-) = ye (22° (50) +8°(50)) 


T(:) “lst (t,)) 


because the life of the fiber is unity, the failure time, t,, 


may be solved for at the stress level S,. 
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APPENDIX B. PROOF TESTING PROCEDURES 


Note: At least 24 hours prior to conducting the fiber 
loading, power up the HP-85 computer, HP-3497A Data 
Acquisition and the HP-62168 power supply (for the load cell) 
to allow the equipment to stabilize before any data is 
recorded. The excitation voltage to the load cell should be 
set to arrroximately 7.5 volts on the power supply dial, The 


voltage c:n be fined tuned using the HP-3497A as a readout. 


1. Turn on the power supply to the elevator. 

2. Turn the plotter. 

3. Ensure that the drive plunger and the load plunger 
on the elevator operate smoothly and there are no 
air bubbles trapped in the plungers or the 
connecting line. 

4. Ensure that there is a full supply of paper for both 
the HP-85 printer and the plotter. Also use fresh 
pens for the plotter. 

5. Perform load cell calibration. 

Overview: The calibration procedure uses a 
program written by LT Bell to determine the slope 
and intercept of the calibration curve for the 


150 gm load cell. Calibration weights are used 
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I. 


in even steps sizes to obtain a plot of the 

calibration curve. Weights ranging from 5 grams 

to 25 grams are normally used in 5 gram steps. 

Insert Bel] Fiber Test cassette in HP-85. 

Type clear - this clears the screen. 

Type load "LDCALB" - the computer searches the tape 

for the program. 

Type run 

Answer the prompts on the screen. 

Enter number of calibration points required (usually 

5). 

Place calibration load on center of load cell. 

i Enter load level. 

ii. When system stops reading data remove weight. 

iii. Repeat step G for each calibration point. 

Enter plot axes data. 

p ae Enter maximum load to be plotted on calibration 
curve. 

ii. Enter minimum load to be plotted on calibration 
curve. 

iii. Enter maximum x value - (read off from tape-- 
approx. -1.0). 

iv. Enter minimum x value (read off from tape-- 
approx. -5.0). 

If large plot is desired, set up plotter and follow 


cues on screen. 
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J. Repeat from step 5.D several times and compare the 
calibration coefficients obtained. The slope coefficient 


A will be used in the program LOADS. 
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A. 


10 
20 
30 
40 
50 
60 
70 
80 
90 
100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 


APPENDIX C. PROOF TESTING SOFTWARE 


LOADCELL CALIBRATION PROGRAM 


! ** LDCALB** 
! JIM NAGEOTEE OCT 26, 1986 
CLEAR 
DIM A(100), B(100), Z1(120), Z2(120) 
DISP *"LOAD CELL CALIBRATION" 
C’ BAR 709 
DISP "NUMBER OF CALIBRATION POINTS, Nil=" 
INPUT Ni 
FOR K = 1 TO Ni 
DISP "INPUT LOAD LEVEL" 
INPUT L(K) 
W = 23 @ E=#5000 
REM :OUTPUT 709; "ARVRIVT3VDS5VAOAE2" 
ON TIMER# 1, E GOTO 240 
Iel 
OUTPUT 709; "AIO" 
ENTER 709; A(I) 
WAIT 10 
OUTPUT 709; "AI10" 
ENTER 709; B(I) 
IeztTI +1 
WAIT 23 
GOTO 160 


240 OFF TIMER # 1 


250 
260 


NeI-1 
Pl=0 @ P2=0 @ Qi=0 @ Q2=0 


270 X1=0 @ X2=0 @ BO = 0 


280 
290 


FOR I=N 
BO=BO+B (TI) 


300 X1=#=X1+A(I) /B(TI) 


310 
320 


X2=#2X2+(A(I)/B(I))*2 
Pl=#P1+A(I) 


330 P2=P2+A(I)%*2 
340 Q1=Q1+B (I) 
350 Q2=Q2+B(I)“*2 


360 
370 
380 


NEXT I 
E9=P1/N 
D9=P2/N-E9“*2 


390 S9=mE9*.5 
400 C9=#ABS (S9/E9*100) 


410 
420 





E8=Q1/N 
D8=Q2/N-E8*2 
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430 
440 
450 
460 
470 
480 
490 
500 
510 
$20 
530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 
770 
780 
790 
800 
810 
820 
830 
840 
850 
860 
870 
880 
890 
900 
910 
920 
930 





S8=D8* .5 

C8=ABS (S8/E8*100) 

X5 (K) =X1/N 

BO=B0/N 

V(K) #X5 (K) 

X6=2X2/N 

D=X6-X5 (K) *2 

S1(K) =D*.5 
C3=£S1 (K) /X5 (K) 

C3=C3*100 

C3=£ABS (C3) 

PRINT "LOAD =’ ;L(K) ; "gm" 
PRINT "MEAN We";X5 (K) ; "V/V" 
PRINT "S.D. V=";S1(K) ; "v/v" 
PRINT "C.V. (&)=";C3 
PRINT "L.C. OUTPUT VOLTAGE " 
PRINT "MEAN V=";E9; "v" 
PRINT "S.D. V=";S9: "vy" 
PRINT "C.V. (%)=";C9 
PRINT "EXCITATION VOLTAGE " 
PRINT "MEAN V=";E8; "Vv" 
PRINT "S.D. Ve";S8: "Vv" 
PRINT "C.V. (%)=";C8 
PRINT "" 

NEXT K 

X1=0 @ Y1=0 

X3=0 @ X2=0 

FOR K=1 TO N 

X1=“X1+V(K) 

Y1leY1=+L (K) 
X3=X3+V (K) *L (K) 
X2~X2+V(K) “2 

NEXT K 

X1=X1/N1 

X2=X2/N1 

X3=2X3/N1 

Yl=#Y1/N1 

D1=#X2-X1“*2 
Cl=(X3-X1*Y1) /D1 
C2=(X2*Y1-X3*X1) /D1 
PRINT "L=A*V+B" 

PRINT "A= ";C1 

PRINT "B= ";C2 

R6=0 

FOR K=1 TO Ni 

G2= (L(K) -C1*v(K) -C2)“*2 
R6=R6+G2 

NEXT K 

PRINT "RESIDUAL = ";R6 
DISP "SCREEN GRAPH ROUTINE" 
DISP "MAX. LOAD ="; 
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940 INPUT Y2 

950 DISP "MIN LOAD ="; 
960 INPUT Y1 

970 U2=#X5 (1) 

980 UleX5 (1) 

990 FOR I=1 TO Nl 


1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1150 
1160 
1170 
1180 
1190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1440 


IF X5(I)>U2 THEN U2=X5 (I) 
IF X5(I)<U1 THEN U1=X5 (I) 
NEXT I 

PRINT *"X MAX (1000V/V)= ";U2*1000 
PRINT "X MIN (1000V/V)= ";U1*1000 
DISP "INPUT MAX FOR X= "; 
INPUT U2 

DISP "INPUT MIN FOR X= "; 
INPUT U1 

FOR I=#1 TO Nl 

B1 (I) #=X5 (I) #1000 

B2 (I) #=S1 (I) #1000 

NEXT I 

S1#160/ (U2-U1) 

S2=140/ (Y2-Y1) 

FOR I#1 TO 5 

GRAPH @ CLEAR 

LDIR 0 

SCALE -48,208,-36,156 
XAXIS 0,32,0,160 

YAXIS 0,28,0,140 

PEN 1 

PENUP 

W1i=160/ (U2-U1) 

W2#140/ (Y2-Y1) 

FOR I=0 TO 5 

U=xU1+ (U2-U1) /5*I 
P3=(U-U1) *W1 

MOVE P3,-10 

LABEL VALS (U) 

NEXT I 

FOR I#0 TO 5 
YxY1+(¥2-Y1) /5*I 
Q3=(Y-Y1) *W2 

MOVE -23,Q3 

LABEL VALS (Y) 

NEXT I 

MOVE 40, -23 

LABEL "1000*Vout/Vexc" 
DEG 

MOVE -29,25 

LDIR 90 

LABEL "Force (gm) " 

LDIR 0 

FOR I#1 TO Ni 
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1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 


1900 
1910 
1920 
1930 
1940 


21 (I) = (B1 (I) -U1) *wi 

22 (1) =(L(I) -¥1) *W2-2 

MOVE 21(I) ,22(T) 

LABEL *** 

NEXT I 

Gi=#C1/1000 

FOR I=0 TO 101 

UsU1+ (U2-U1) /100*I 

Y¥=G1*U+C2 

IF Y<Y1 OR Y>Y2 THEN 1580 

Z1 (I) =(U-U1) *W1 

Z2 (I) = (¥-Y1) *W2 

PLOT Z1(I) ,Z2 (I) 

NEXT I 

COPY 

DISP "DO YOU WANT A HARD COPY ON PLOTTER? (Y/N) " 
INPUT K8$ 

IF K8S$="N" THEN 2170 

PRINT IS 10 

CONTROL 10,5 ;48 

OUTPUT 10 ;"IN;SP1;IP2400,1,600,8800,6900,: 
OUTPUT 10 ;"SC0,1000 ,0,1000" 
OUTPUT 10 ;"PU0,0PD0,1000,1000,1000,1000,0,0,0PU" 
W3=1000/ (U2-U1) 

W4=21000/ (Y2-Y1) 

OUTPUT 10 ;"SI0.2,0.3;TL1.5,0" 

FOR I=#0 TO 5 

YeY1+ (Y2-Y1) /5*1I 

Y4= (Y-Y1) *wW4 

Y=INT(Y) 

Y4=INT (Y4) 

OUTPUT 10 ;"PA 0,",Y4,"YT;" 

OUTPUT 10 ;"CP-5.-0.07;LB"; Y;CHRS (3) 
NEXT I 

FOR I=#0 TO 5 

UsU1+ (U2-U1) *I/5 

U4= (U-U1) *W3 

U4zn INT (U4) 

UszxINT (U) 

OUTPUT 10 ;"PA";U4,",0;XT;" 

OUTPUT 10 ;"CP-1.3,-1;LB";U; CHRS (3) 
NEXT I 

OUTPUT 10 ;"SI.30, .42" 

OUTPUT 10 ;"PA400,0;CP-2,-2.3;LB1000*Vout/Vexc" ; CHRS (3) 
OUTPUT 10 ;"PA0,460;DI0,1;CP-2.6,2.6;LB FORCE (GM)"; 
CHRS (3) 

OUTPUT 10 ; "DI; PU" 

FOR I =1 TO N1 

Z1(I) =(B1(I) -U1) *W3 

Z2 (I) = (L(I) -¥1) *W4 

Z1 (I) «INT (Z1 (TI) ) 
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1950 22 (I) = INT (22 (TI) ) 


1960 


CHRS (3) 


1970 NEXT I 
1980 OUTPUT 10 :*PU" 

1990 FOR I=0 TO 101 

2000 UsU1+ (U2-U1) /100*I 

2010 Y=G1*U+C2 

2020 IF Y<Y1 OR Y>Y2 THEN 2060 

2030 21(I)=(U-U1) *W3 

2040 22(I)=(Y-Y1) *wa 

2050 OUTPUT 10 ;"PA";Z1(I),Z2(I):"PD" 
2060 NEXT I 

2070 OUTPUT 10 ;"PU0,900,100,900" 
2080 OUTPUT 10 ;"PU*" 

2090 OUTPUT 10 ;"SI.22,.38" 

2100 DISP "ENTER THE LEGEND. ENTER ‘0’ 
2110 INPUT P7$ 

2120 IF P7$="0" THEN 2150 

2130 OUTPUT 10 ;"CP;LB";P7$;CHRS$ (3) 
2140 GOTO 2100 

2150 PRINTER IS 2 

2160 DISP "END LDCALB" 


2170 END 
PROOF TEST PROGRAM 
10 ! 
15 ! INCORPORATES PLOT ROUTINE* 
20 ! JIM NAGEOTTE SEP 6,91 
30 ! 
40 CLEAR @ DISP @ DISP 


45 


50 
55 
56 
57 
60 
70 
80 


85 


90 
120 
130 
140 
150 


DISP "FOR USE ON THE 150 GM 

CELL! * 

DISP " QUIKLD PROGRM" 

A=-12621 

B=-23 .291198677 

BEEP 

PRINTER IS 2 

SHORT A(55),B(55) ,P(600) 

PRINT "IF PROGRAM HALTS TYPE 
‘CONT 100’ TO RETURN TO MEN 

u*" 

PRINT @ PRINT “MAKE SURE THA 

T LINE #251 IS CURRENT." 

PRINT @ PRINT 

ON KEY# 2,*ADJ B" GOSUB 1000 

ON KEY# 1,"DATE " GOSUB 210 

ON KEY# 3, "LOAD" GOSUB 1320 
! ON KEY# 3, “WEIGH" GOSUB 10 

90 
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OUTPUT 10; "PAO,460;DI0,1;CP-2.6;2.6;LB FORCE (gm) 


TO EXIT" 











151 ! ON KEY# 5S, "MX LOAD" GOSUB 
4000 

155 ON KEY# 4, "PLOT" GOSUB 3000 

160 ! ON KEY# 8, "END" GOTO 190 

170 KEY LABEL 

180 GOTO 120 

190 CLEAR @ DISP @ DISP @ BEEP @ 
DISP * END OF HPLOR" 

200 END 

210 ! * INPUTS * 

220 CLEAR @ DISP “ENTER THE DATE 
a 


230 INPUT DS 

251 Wa4.8 

280 CLEAR 

290 RETURN 

1000 ! * ADJUST B FOR 0 * 

1010 CLEAR @ DISP @ DISP 

1020 DISP "APPLY ZERO LOAD TO LO 
AD CELL" 

1030 BEEP @ DISP 

1040 GOSUB 1130 

1050 Bl=#P1/23 

1060 BeB-Bi 

1070 GOSUB 1185 

1080 CLEAR @ RETURN 

1090 ! * WEIGH * 

1100 ! SHORT ROUTINE TO CHECK 

1110 ! ACCURACY OF CALIBRATION 

1120 CLEAR 709 @ CLEAR 

1130 OUTPUT 709 ;"AI10" 

1140 ENTER 709 ; V 

1150 Al=A/V 

1170 BEEP 

1185 Pi=0 

1190 FOR I=1 TO 25 

1200 OUTPUT 709 ; "AIO" 

1210 ENTER 709 ; L 

1220 PeL*Al+B 

1230 PINT (P*1000) /1000 

1235 IF I <3 THEN 1260 

1240 DISP P, 

1250 Pil=#=P1+P 

1260 NEXT I 

1270 DISP 

1280 PRINT "Vout=";V 

1290 PRINT "LOAD ="; INT(P1/23*10 
00) /1000;" gm" 

1300 BEEP 

1310 RETURN 

1320 ! * LOAD IT * 
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1390 


1400 
1403 
1405 
1407 
1409 
1410 


1420 


1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1542 
1543 
1544 


1550 
1560 


1570 
1580 
1585 


1590 
1600 
1610 
1620 
1630 
1640 


1645 
1650 
1660 


1670 
1680 
1690 
1700 
1750 
1760 


1770 


CLEAR @ PRINT @ PRINT @PRI 
NT @ PRINT 

PRINT DS 

DISP "ENTER THE SAMPLE#" 
INPUT S2$ 

DISP "ENTER THE STATION #" 
INPUT S1 

DISP @ DISP @ DISP “PREPARE 
LOADER" 

DISP @ DISP "HIT ‘CONT’ WHE 
N READY." 

PAUSE 

BEEP @ CLEAR 

CLEAR 709 

OUTPUT 709 ; AI10" 

ENTER 709 ; V 

Al=A/V 

OUTPUT 709 ; "AIO" 

ENTER 709 ; L 

P=#«L*A1+B 

P#ABS (P*#100/100) 

Ps INT (P*1000) /1000 

BEEP 

PRINT @ PRINT "WEIGH MEASU 
RES"; P 

DISP "WEIGH IS ";P 

DISP "START LOADER AND HIT 
‘CONT’ SIMULTANEOUSLY" 
PAUSE 

C9=0 @ I=0 @ P3=W 

DISP @ DISP @ DISP " LOAD 
ING" 

IeI+1 

OUTPUT 709 ; "AIO" 

ENTER 709 ; L 

P(I) =L*A1+B 

IF I<4 THEN 1590 

! IF P(I)>P(I-1) THEN P4=P( 
I) 414.17 

IF P(I)<P3 THEN P3=P(I) 

IF P(I)<W*.9 THEN C9=#1 

! IF C9=x21 AND P(I) >W*.95 TH 
EN GOTO 1690 

WAIT 25 

GOTO 1590 

CLEAR @ DISP @ DISP @ DISP 
BEEP 100,100 
P2=INT ( (W- P3) #100) /100 
DISP "FIBER BROKEN AT ";P2; 


PRINT S2$;" #";S1;" BROKEN 
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1775 
1780 
1840 
1850 
1860 


1863 


1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 
1950 
1960 


1970 
1980 


1990 
2000 
2010 


2020 
2030 


2040 
2045 


2046 


2050 
2060 
3000 
3001 
3003 


3004 
3006 
3010 
3020 
3030 
3040 





AT ";P2; "gms." 

PRINT P2; "gms" 

CLEAR @ RETURN 

! * TIMER * 

T1=TIME 

BBEP @ CLEAR @ DISP "TURN O 
FF LOADER" @ BEEP 

DISP @ DISP @ DISP "TIMING 
FOR TWO MINUTES. * 

OUTPUT 709 ; "AIO" 

ENTER 709 ; L 

PeL*A1+B 

T2=TIME 

T3=0T2-T1 

IF T3>120 THEN GOTO 2000 
IF PoW*.9 THEN GOTO 1950 
GOTO 1870 

! * IT BROKE * 

CLEAR @ DISP @ DISP @ BEEP 
100,250 

DISP "FIBER BROKEN AT ";T3; 
" SEC." 

PRINT "FIBER #";S1;" BROKEN 
AT ";T3;" SECONDS." 

CLEAR @ RETURN 

! * END OF TIME * 


CLEAR @ DISP @ DISP @ DISP 
@ DISP @ DISP @ DISP "FIBER 
OK" 

BEEP 50,200 


PRINT "FIBER #";S1;" INTACT 


PRINT @ PRINT 

BEEP @ DISP 2 DISP "MAKE PL 
OT!" @ BEEP 

DISP @ DISP "HIT CONT" @ PA 
USE 

CLEAR @ RETURN 

END 

! **PLOTTER** 

CLEAR 

BEEP @ FOR T#=1 TO 100 @ NEX 
T @ BEEP @ DISP @ DISP @ 
DISP "INSTALL PAPER IN PLOT 
TER!" 

DISP @ DISP "HIT CONT" 
PAUSE 

CLEAR 

Y2=s5 

X2=I 

NeI 
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3050 
3060 
3070 
3080 
3090 


3100 
3110 
3120 


3130 
3140 
3150 
3160 
3170 


3180 
3190 
3200 
3210 
3220 
3230 


3240 
3250 
3260 


3270 


3280 
3290 
3300 
3310 
3315 
3320 
3330 
3340 
3350 


3360 
3370 
3380 
3390 
3400 
3410 
3420 
3440 
3460 


$3=1000/X2 

$4=1000/Y2 

PRINTER IS 10 

CONTROL 10,5 ; 48 

OUTPUT 10 ;*IN;SP1;1I1P2400,1 
600, 8800, 6900; ”* 

OUTPUT 10 ;*SC0,1000,0,1000 
* 


OUTPUT 10 ;*PU0O.0OPD0,1000,1 
000,1000,1000.0.0.0PU" 
OUTPUT 10 ;"SI0.2,0.3;TL1.5 
,0" 

FOR I=#0 TO 5 

Ys INT (Y2/5*I) 

Y4=INT (Y*S4) 

OUTPUT 10 "PA 0,",Y4, "YT; " 
OUTPUT 10 ;"CP-5,-0.07;LB"; 
Y¥; CHRS (3) 

NEXT I 

FOR I=#1 TO 5 

X=eX2/5*1I 

X42INT (X*S3) 

OUTPUT 10 ;"PA";X4,",0;XT;" 
OUTPUT 10 ;"CP-1.3-1;LB";X 
; CHRS (3) 

NEXT I 

OUTPUT 10 ;"SI.30, .42" 
OUTPUT 10 ;"PA220,0;CP4,-2. 
3;LBData Point"; CHRS (3) 
OUTPUT 10 ;"PA0,460;DI0,1;C 
P-2.6,2.60;LB LOAD (gm) ";CH 
R$ (3) 

PRINT "SPO" 

PRINT "SP2" 

FOR I=1 TO N 

ZisINT (I*S3) 

P52£W-P (I) 

Z2=INT (P5*S4) 

OUTPUT 10 ; "PA", Z1,2Z2; "PD" 
NEXT I 

OUTPUT 10 ;"PpU0,900,100,900 
f 


OUTPUT 10 ; "PU" 

OUTPUT 10 ;"SI.22, .38" 
PRINT *SPO" 

OUTPUT 10 ;"SI.22, .38" 
PRINT "IN;SP1; PA3000,6500" 
DISP "ENTER ONE LINE LABEL" 
INPUT BS 

C$=D$ 

AS$=S2$ 
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3470 
3480 
3490 
3495 
3500 
3510 
3520 
4000 
4005 


4010 


OUTPUT 10 ;"CP;LB";AS;" ON 
#";S1;CHRS (3) 

OUTPUT 10 ; "CP; LB"; BS; CHRS ( 
3) 

OUTPUT 10 ; "CP;LB";CS;CHRS ( 
3) 

! OUTPUT 10 ;"CP;LB"; MAX 
LOAD IS ";P4 

PRINTER IS 2 

CLEAR 

RETURN 

P4eP4+14.17 

DISP "MAX LOAD ON FIBER = " 
3P4 

RETURN 
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APPENDIX D. REDUCED TIME PARAMETERS 
It is assumed that the failure times of the fibers can be modeled 


by a Weibull distribution as 


F(t) =1-exp (- (©) **) 
B. 


where t is the independent variable, 8, is the location parameter 
and a, is the shape parameter. 

To understand the parameters associated with the reduced time, 
consider a stress rupture history where the fiber is imagined to be 
instantaneously loaded to a sustained stress level S, until failure 
at time ty. Using a breakdown rule of the power law form, the 


reduced time is 


*f og 
er | —1)P 
T(+) gf (qorae 
After integration, 
t.,S 
TC) (8 


The reduced time may be substituted into the failure distribution 


as stated in section II.B.3, yielding 
‘ete ~ Cr S,,°\* 
F(+) =1 expl-{(4F) (=) )} 


and set equal to the Weibull distribution above. 
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by comparing exponents, a = a,, and rearranging gives 


Be, S51, °. 
= (5)? 1 


If a deterministic approach is taken, 8, may be replaced by simply 
ts, the failure time. If € is chosen to equal te, then A is the 
intrinsic strength at time €. The parameter p is simply the slope 


of the breakdown function on a log-log plot. 
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APPENDIX E. PROOF TEST STRESS HISTORY 


The proof test history is shown below and is divided into six 


regions as described. 


stress 





Figure 22. Proof Test Stress History 


S(t) = Lyt ;O st «<ty Region I 
= Si ity st < ty Region II 
= Lo (t-t2) + S, ; tg st <ty, Region III 
= So i tz rr < t, Region IV 
= L3 (t-t,) + S2 ; tgs t <te Region V 
=S, ; t >te. Region VI 


Using a breakdown rule of the power form and the stress history as 


above, the fractional life used in each region is the sum of the 
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t. ts 
T()=4 f(A yeaged f(y rae 
[*) t; 


ty t ts ty 
1¢,%,(t-t,),, 1, S2 17, £,(t-t,) 1,51). 
es ea a ee er a ee 9 a) ede 


integrals up to the failure time. 


Performing and evaluating at the limits of integration yields, 





gt (ta A Say? (pc 
T(:) ='> (pziy E a) (t,-t,) 
ta A Ea ey ay) 
€ L,(p+1) A A 
Ss. e 
+4 (2) (t-t,) 
res A (L, (€,-,) +S.) gases ae 
EL, (p+1) A A 


1,51," 
| =) (t,-ts) 
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A. SIMDATA 


% PROGRAM SIMDATA 

Written by LT Greg Morin 

% This program produces simulated data which can be used for 
analysis of 

% of a Weibull distribution model. Inputs include the desired 
population 

% size and the underlying shape and location parameters, N, 
alpha, and beta 

% respectively. Outputs include the following column vectors: 


% X= [population] 

% xi=[(exact data] 

% xr=[right censored data] 

% Fstarx=(underlying F* values for population] 
+ Fstar=(Expected Ranking F* values] 

clear 

clg 

Nexinput (’Enter the population size: ‘); 


a=input (‘Enter desired underlying alpha: ")3 
beinput (‘Enter desired underlying beta: "); 

% Function to simulate data and plot fstar for expected and true 
rank 

% of the population. 

[X, Fstar, Fstarx] =populat (N,a,b) ; 

hold on 

% Function to determine exact and right censor data 
(xi, Fstrc,m,xr] =rightcns (X,N,Fstar) ; 

hold off 

save xi xi /ascii 

save xr xr /ascii 

save X X /ascii 

save Fstarx /ascii 

save Fstar /ascii 


B. FRACTIONAL LIFE 


Eteteeteee Program Name: Eqvl_1, written for MATLAB ****#ttex 
Written by LT Joe Woodward 


clear %CLEAR ALL VARIABLES IN WORKSPACE 
format long e &%&SET OUTPUT FORMAT TO 16 DIGITS 
% VARIABLES DEFINED 
% 
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16 


PC PHRMA AHHHAHAHA PAH H EHH ARH HEP BHM 
E 
i) 
m 





intrinsic strength of the fiber 
stress level that fiber was realized at during 
ramp loading 
the size of the sample 
the initial loading rate used on the fibers 
(gm/sec) 
the loading rate used on the proof test (gm/sec) 
the unloading rate used on the proof test (gm/sec) 
intrinsic time used to determine strength 
the slope of the breakdown law 
sustained stress level 
stress level of proof test 
time sustained stress level is reached 
time proof test begins 
time proof load is attained 
time proof load is removed 
time sustained stress level is reached after 
proof load 
duration of proof load 
failure time of non-proof tested fibers 
failure time of proof tested fibers 
non-proof tested fiber failure stress 
proof tested fiber failure stress 


arrays that store the fractional life consumed 
during each region of the proof test stress 
history. 1 is how much life was used during the 
initial loading, 11 the life used up to proof 
loading at the sustained stress level, etc. 

the total life consumed by each fiber. Will equal 
rnity when the fiber fails. 


Load the simulated realized stress values 


~ 
oO 
i) 
a 
* 


% The following section asks for inputs from the operator for use 


* in the 
ldoti = 


idot2 = 
idot3 = 
rho = 


teehat = 
$1 = 


program. 

input (’Enter the initial loading rate for the fibers 
(gm/sec), ldoti: ‘); 

input (‘Enter the loading rate for the proof test 
(gm/sec), ldot2: '); % 

input (‘Enter the unloading rate for the proof test 
(gm/sec), ldot3: ‘); 

input (‘Enter the slope of the breakdown law, rho: '); 
input (‘Enter the intrinsic time, teehat: ‘); 

input (‘Enter desired sustained stress level, Si: '); 


% CALCULATE THE TIME THAT THE SUSTAINED STRESS IS REACHED 
timel = Si/ldoti; 
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% DETERMINE THE SAMPLE SIZE (NUMBER OF X VALUES LOADED) 
N=length (x) ; 


% THIS LOOP FINDS THE INTRINSIC STRENGTH OF EACH SAMPLE, THE LOOP 
% IS REPEATED WN NUMBER OF TIMES 
for ie1:N 

A(i) = ((x(i)* (rho+1)) / (teehat*ldot1* (rho+1)))*(1/rho) ; 


% FIND THE FAILURE TIME OF EACH FIBER DURING RAMP LOADING 
t_ramp(i) = x(i)/ldoti; 
end 


% SET THE INTRINSIC LIFE OF EACH FIBER TO UNITY 
tauel1.0; 





% CLEAR ANY GRAPHS AND THEN PLOT THE RAMP HISTORY (FAILURE TIME % 
% VS. LOAD) 

clg 

plot (t_ramp,x,'rX’) ,grid 

xlabel (‘time (sec)’),ylabel(’load (gm) ’) 

title(’Fiber Failure Under a Constant Loading Rate’); 

pause 


% THIS LOOP IS FOR THE NON-PROOF TESTED HISTORY 
for j = 1:N 
tnptf = t_ramp(j); 
sf = ldoti*tnptf; 


IF THE FAILURE TIME OF THE FIBER IS GREATER THAN THE TIME AT 
WHICH SUSTAINED STRESS LEVEL IS REACHED, THEN COMPUTE THE 
INTRINSIC STRENGTH AND THE FAILUR TIME GIVEN THAT THE FIBER 
SURVIVED AT LEAST UNTIL THE SUSTAINED STRESS LEVEL. THE STRESS 
LEVEL AT WHICH THE FIB. FAILED AT IS THE SUSTAINED STRESS 
LEVEL, 81. . 
if tnptf>=«timel 
a=(((1dot1/A(j) ) *rho) * ((time1* (rho+1) ) / (rho+1) )); 
tnptf = (tau*teehat-a)*((A(j)/S1)*rho) + timel; 
sf = S1; 
end 


OHHH 


% FILL THE ARRAYS TNPT(J) AND S(J) WITH THE FAILURE TIMES AND 
% CORRESPONDING STRESSES FOR PLOTTING. 

tnpt (j) «tnptf; 

s(j) = sf; 

end 


% PLOT THE NON-PROOF TESTED FIBERS ON SEMILOG (IN TIME) VS. 
% FAILURE STRESS 

semilogx(tnpt,s,’y+’ ) 

xlabel(’log time (sec)’),ylabel(’stress’),grid 

hold 
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pause 


% THIS SECTION ASKS FOR INPUT REGARDING THE PROOF TEST. 

as = input (‘Enter proof load magnitude, Ds: '); 

time2 = input(’Enter the time the proof load is conducted (sec), 
time2: '); 

time6é = input(’Enter the duration of the proof load (sec), 
time6: '); 


% COMPUTE THE TIME THE PROOF LOAD IS REACHED (TIMES). 
time3 = time2 + dS/ldot2; 


% COMPUTE THE TIME THAT THE UNLOADING OF THE PROOF TEST STARTS 

% (TIME4), WHERE TIMEG IS THE DWELL TIME, i.e., THE DURATION THAT 
% THE PROOF LOAD IS MAINTAINED. 

time4 = time3 + time6; 


% COMPUTE THE TIME THAT THE ORIGINAL SUSTAINED STRESS LEVEL IS 
% REACHED AFTER UNLOADING THE PROOF LOAD (TIMES). 

% proof load 

timeS = time4 + (Ds/(-ldot3)); 


% S2 IS THE MAGNITUDE OF THE PROOF LOAD, ORIGINAL SUSTAINED 
% STRESS LEVEL PLUS THE INCREMENTAL PROOF LOAD. 
S2 = S1 + ds; 


% THIS LOOP IS FOR THE PROOF TESTED FIBERS. 
for j = 1:N 

tptf = t_ramp(j); 

ssf = ldotittptf; 


Ee eeeRetearetetaneaeeRREGION TESS SAeAAAKeAEKAERaAAKEEHaARUHaaakEne 


% IF THE FAILURE TIME OF THE FIBER IS LESS THAN THE TIME AT 
% WHICH THE SUSTAINED STRESS LEVEL IS REACHED, THEN COMPUTE THE 
% FRACTIONAL LIFE, (1(4j)), CONSUMED GIVEN THAT THE FIBER FAILED 
% ON LOADING. 
if tptf<=timel 
1(j)=(1/ (teehat* (rho+1) )) *((1ldot1/A(j))“*rho)... 
* (tptf* (rho+1)); 
else ; 


% IF THE FAILURE TIME IS GREATER THAN THE THE TIME AT WHICH THE 
% SUSTAINED STRESS LEVEL IS REACHED, THEN COMPUTE THE FRACTIONAL 
% LIPE BASED ON TIME] 
1(j) =(1/ (teehat* (rho+1) )) * ((1dot1/A(j))*rho)... 
* (time1* (rho+1) ); 
end 
16 (j) #1(4); 
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Beeeoeaseeeaeeeeee PEGI ON [Ie eseeeeeterseseesseaesesaeeesaaeesares 


%& IF THE PAILURE TIME IS GREATER THAN THE TIME AT WHICH THE 
& SUSTAINED STRESS LEVEL IS REACHED, THEN COMPUTE THE FAILURE 
% TIME GIVEN A NEW STRESS HISTORY, i.e., THE FIBER SURVIVED 
% LOADING AND SPENT SOME OF ITS LIFE UNDER A SUSTAINED CONSTANT 
% LOAD, Sl. THE FAILURE TIME IS COMPUTED BY SUBTRACTING THE 
% FRACTIONAL LIFE CONSUMED DURING THE INITIAL LOADING PROCESS 
& PROM THE TOTAL LIFE TAU. 
if tptfotimel 
tptf = (tau-1(4j))*(teehat*(A(j)/S1)“rho) + timel; 
ssf = Sl; 


if tptf<=stime2 
11 (5) =(1/teehat) * ((S1/A(j))“*rho) * (tpt£-timel) ; 
else 
11(j)=(1/teehat) * ((S1/A(j))“*rho) * (time2-time1) ; 
end 
16 (j) #1 (5) +11(5); 


RACHA HAEAEAECHHERREGCION [IOC eee wee ee aeeaeaeeeeeeeeaaeeaeaannen 


EACH IF STATEMENT CHECKS TO SEE IF THE FAILURE TIME FALLS 
BETWEEN THE START OF A NEW REGION AND THE START OF THE NEXT 
REGION. IN THIS CASE, IS THE FAILURE TIME BETWEEN THE START OF 
THE SUSTAINED STRESS REGION AND THE START OF THE PROOF LOADING 
REGION. IF IT IS, THE FAILURE TIME IS COMPUTED BASED ON THE 
GERMAIN STRESS HISTORY AND THE FRACTIONAL LIFE IS COMPUTED 
BASED ON THAT FAILURE TIME. IF THE FAILURE TIME IS GREATER 
THAN WHEN THE NEXT REGION STARTS (IN THIS CASE THE PROOF 
TESTING) THAN THE FRACTIONAL LIFE CONSUMED IS BASED ON THE 
ABSOLUTE TIME SPENT UNDER THE SUSTAINED LOAD. 
if tptf>time2 
tptf=(((tau-1(j)-11(j)) * (teehat*ldot2* (rho+1) /A(j))... 
+((S1/A(j))* (rho+1)))*(1/ (rho+1)) - (S1/A(j)))... 
*(A(j) /ldot2) +time2; 


HHHHLHLHHE eS 





% THE FAILURE STRESS WILL BE BETWEEN THE ORIGINAL SUSTAINED 
% STRESS LEVEL, S1 AND THE PROOF TEST LEVEL S2. 
ssf = ldot2*(tptf-time2)+Si;'. 


if tptf<=time3 
12 (j) =(1/teehat) * ((((1ldot2* (tptf-time2)... 
+S1) /A(j))* (rho+1) - ((S1/A(j))* (rho+1))) *(A(5)/... 
(ldot2* (rho+1)))); 
else 
12 (j) =(1/teehat) * (( ( (ldot2* (time3-time2)... 
+81) /A(j))*(rho+1)... 
- ((S1/A(5))* (rho+1))) *(A(j) / (ldot2* (rho+1)))); 
en 
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16 (j) #1 (3) #115) +12 (5); 


& SHAaeetaneananantee REGION LVEF eeeeeeeeeeeeeeaeaeeaeaaeaeeeceaeaarenn 


if tptfo-time3 
tptf =(tau-1(j)-11(j)-12(j))* (teehnat*(A(j)... 
/S2)“rho) +time3; 
ssf = S2; 


if tptf<=time4 
13 (j)=(1/teehat) * (((S2/A(j)) “*rho) * (tpt£-time3) ) ; 
else 
13 (j) =(1/teehat) *(((S2/A(j)) “rho) * (time4-time3) ) ; 
a 


en 
16 (5) =1 (5) #11 (5) +12 (5) +13 (5); 


ReteneceeeranaeereaRRGION Vere eaeseaeeeseaeaeaeeeeeeeaeeeeacene 


if tptf>-time4 
tptfe- (( (tau-1(j)-11(j)-12(5)-13(5))... 
* (teehat* (-1dot3) * (rho+1) /A(j))+((S2/A(j))... 
“(rho+1)))*(1/ (rho+1) ) - (S2/A(j))) *(A(5)/... 
ldot3) +time4; 
ssf = S2+ldot3* (tptf-time4) ; 


if tptf<\etime5 
14 (j) =- (1/teehat) * ((( (ldot3* (time4-tptf)... 
+S2) /A(j))* (rho+1) ) - ((S2/A(j))*(rho+1... 
)))*(A(3j) / (1dot3* (rho+1))); 
else 
14(j) =- (1/teehat) * (( ( (ldot3* (time4-timeS)... 
+S2) /A(j))* (rho+1)) - ((S2/A(j))*(rho+1... 
)))*(A(3) / (ldot3* (rho+1))); 
end 
16 (5) =1 (3) +11 (5) +12 (5) +13 (5) 414(5); 


RBACeHeaeaeeaanareaaaeaeaREGION VI St eeseAeaeeaneaaaaaeaaekekaadeakaneanank 


if tptf>time5 
tptf=(tau-1(j)-11(j) -12(5) -13(j)-14(5))... 
*(teehat* (A(j)/S1)“*rho) +timeS; 
ssf= Si; 
15(j)=(1/teehat) *((S1/A(j)... 
)*rho) *(tptf-timeS) ; 
16 (5) #1 (5) +11 (5) 412 (5) +13 (5)... 
+414 (j)+15(j); 
end 
end 
end 
end 
end 
tpt (j)=tptft; 
ss (j)=ssf; 
end 
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% PLOT THE PROOF TESTED STRESS HISTORY, FAILURE TIME VS. STRESS 
semilogx(tpt,ss,’rxX’) 

save SS.m ss -ascii 

save TPT.m tpt -ascii 

save S.ms -ascii 

save TNPT.m tnpt -ascii 

format short e 


% THE MATRIX T IS THE FRACTIONAL LIFE OF EACH FIBER CONSUMED 
% DURING EACH REGION AND THE FAILURE TIME. 

Tx(1’ 12’ 12’ 13’ 14’ 15’ 16° tpt]; 

save T.txt T -ascii 

save T.m T 


Cc. PARAMETERS USED IN DETERMINISTIC LIFE STUDY 


For all runs, a = 5, B = 20, p = 40, € = 1 and tz = 50000. 


Run #1: L,= .8, Lo= .8, L3= -.8, S; = 14.1, Sz = 18.6, 
taweil =1- 


Run #2: Ly= .8, Lo= 50, L3= -50, S, = 14.1, S = 18.6, 
Caweil =1- 
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APPENDIX G. 


A. POST-PROOF TEST SOFTWARE 

% eeaesexene Post-Proof Test Software *tt*tttet 

Written by LT Joe Woodward 

clear 

format long e 

% VARIABLES DEFINED 

% 

% 

* x stress level that fiber was realized at during 
% ramp loading 

on the size of the sample 

& Ildot the loading rate used on the fibers (gm/sec) 
* teehat intrinsic time used to determine strength 

% rho the slope of the breakdown law 

* si sustained stress level 

* S2 stress level of proof test 

%* timel time sustained stress level is reached 

* time2 time proof test begins 

* time3 time proof load is attained 

* time4 time proof load is removed 

* time5 time sustained stress level is reached after 
t proof load 

* time6é duration of proof load 

n = 2000;%input (‘Enter the number of samples to be 
calculated, N: ‘); 

alpha = 5;%input(’Enter desired underlying alpha: oe 
beta = 20;%input (’'Enter desired underlying beta: '); 
ldoti «= .8;%input(‘Enter the initial loading rate for the fibers 
(gm/sec), ldoti: '); 


ldot2 = .8;%input(’Enter the loading rate for the proof test 
(gm/sec), ldot2: '); 

ldot3 = -.8;tinput(’Enter the unloading rate for the proof test 
(gm/sec), ldot3: '‘); rs 

rho = 40;tinput (’Enter the slope of the breakdown law, rho: 
")4 


teehat = 1;%input(’Enter the intrinsic time, teehat: ‘); 

$1 = 7;tinput(’Enter desired sustained stress level, S1: '); 
Ds = 3;%input(’Enter proof load magnitude, Ds: '); 

time2 = 50000;%input(’Enter the time the proof load is conducted 


(sec), time2: ‘); 
time6 = 1;%tinput(’Enter the duration of the proof load (sec), 
time6: '); 
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% time that sustained stress level is reached 
timel = S1/ldoti; 


% time that proof load is reached 
time3 = time2 + Ds/ldot2; 


* time that unloading of proof load starts 
times = time3 + timeé; 


% time that original sustained stress level is reached after 
unloading 

* proof load 

timeS = time4 + (Ds/(-ldot3)); 

S2 = $1 + Ds; 

erreetetoutside loop 

deltime=zeros (1,50) ; 

for jj=#1:50 

LL=1; 


% Find the intrinsic strength for each fiber at teehat 
A = ((x*(rho+1))/(teehat*ldoti1* (rho+1)))*(1/rho) ; 
% Find the failure time of each fiber during ramp loading 
t_ramp = x/ldoti; 
* Set the reduced time (Life) of each fiber 
tau = 1.0; | 
% *##** NON-PROOF TESTED FIBERS **** 
tnptf = t_ramp; 
sf = ldotittnptf; 
if tnptf<=timel 
1x=(1/ (teehat* (rho+1) ) ) * ((1dot1/A) “rho) * (tnptf* (rho+1) ); 
else 
lx= (1/ (teehat* (rho+1) ) ) * ((ldot1/A) “rho) * (time1* (rho+1) ); 
eae 


if tnptf>timel 














tnptf = (tau-1x)*(teehat*(A/S1)“rho) + timel; 


sf = Si; 
lix=(1/teehat) * ((S1/A) *rho) * (tnptf-timel) ; 
16x (j) #1x+11x; 

end 

lifex(j)=tnptf; 


% ‘eee PROOF TESTED FIBERS **** 


tptf = t_ramp; 
ssf = ldoti*tptf; 


if tptf<=timel 
L=(1/(teehat* (rho+1) )) * ((1dot1/A) *rho) * (tptf* (rho+1) ) ; 
else 


l=(1/(teehat* (rho+1))) * ((1dot1/A) “rho) * (time1* (rho+1)); 
end 
16(j)=1; 


if tptf>timel 
tptf = (tau-1)*(teehat*(A/S1)“rho) + timel; 
ssf = Sl; 


if tptf<=time2 
11=(1/teehat) * ((S1/A) *rho) * (tptf£-timel) ; 
else 
li=(1/teehat) * ((S1/A) “*rho) * (time2-time1) ; 
end 

16 (3) =1+11; 


if tptf>time2 


tptf=(((tau-1-11) *(teehat*ldot2* (rho+1) /A)... 
+((S1/A) * (rho+1)))*(1/ (rho+1))-(S1/A))... 
* (A/1ldot2) +time2; 

ssf = ldot2* (tptf-time2) +S1; 


7f£ tptf<=time3 


12=(1/teehat) * ( (( (1dot2* (tptf£-time2) +S1) /A)*(rho+1)... 


- ((S1/A) * (rho+1)))* (A/ (ldot2* (rho+1)))); 
else 


12=(1/teehat) * (( ( (1dot2* (time3-time2)+S1) /A)*(rho+1)... 


- ((S1/A) * (rho+1) )) * (A/ (ldot2* (rho+1)))); 
end 
16 (3) =1+11+12; 
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if tptf>-time3 


tptf = (tau-1-11-12) * (teehat* (A/S2) “rho) +time3; 
ssf = S2; 


if tptf<=time4 
eo tereebee Caan) rho) * (tptf£-time3) ); 


«Se tu) ceahac)® CL le3/A\*one)® (eines eames) 


en 
16 (5) #=1411+12+13; 


if tptfotime4 

tptf=- (( (tau-1-11-12-13) * (teehat* (-1dot3) . 
* (rho+1) /A) + ((S2/A) * (rho+1)))* (1/(rho+1)). 
- (S2/A) ) * (A/ldot3) +time4; 

ssf = S2+ldot3* (tptf-time4) ; 


if tptf<=time5 
14=- (1/teehat) * ( ( ( (1dot3* (timed - tptf)+S2) /A)*(rho+1)). 
- ((S2/A) * (rho+1) )) * (A/ (ldot3* (rho+1) )); 


else 


14=- (1/teehat) *((( (1dot3* (time4-timeS)+S2)/A)*(rho+l1))... 
- ((S2/A) * (rho+1) )) * (A/ (1dot3* (rho+1))); 


end 
16 (3) =14+11+12+13+14; 


if tptf>time5 
tptf =(tau-1-11-12-13-14)*(teehat*(A/S1)... 
“rho) +time5; 
ssf= Sil; 
life(j)=tptf; 
15=(1/teehat) * ( (S1/A) * enone GREES time5) ; 
16 (j) =#14+11412+13+14+15; 


end 
end 
end 
end 
if tptf>time5 
ftnew (LL) =tptf; 
ftminemin (ftnew) ; 
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end 


tnpt (j) =tnptf; 
s(j) = sf; 

tpt (j) =tptf; 
ss(j) = ssf; 


end 


tnpt=tnpt’; 
8=s'; 
tpt=tpt’; 
ss=ss’ ; 
Geltime (jj) =(ftmin-time2) ; 
end 
% 
% Routine to develope histogram for no failure time. 
% 


% 

sort (deltime) ; 
numberl=length (deltime) ; 
tmaxl=max (deltime) ; 
tmax1=10910 (tmax1) ; 
tminl=min (deltime) ; 
tmini=#log10 (tmin1) ; 
atl=(tmaxl-tmin1) /10; 
yi=zeros (11,1); 
xtl=zeros (11,1); 


for kel:numberl 
i=1; 
atime=deltime (k) ; 
dtime#1ogi10 (dtime) ; 
for timl=tmini1:dt1l:tmaxi 


tnextletiml+dti; 


if dtime>=timl 
if dtime<tnextl 
y1(i) sy1 (i) +1; 
xt1 (i) #10* ( (timl+tnext1) /2); 
end 
end 
iwi+l1; 
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end 


end 


clg 

semilogx(tnpt,s, ’yX’) 

xlabel(’log time (sec)’),ylabel(’load (gm)’),grid 
title(‘Non-Proof Tested Failure Times’ ) 

pause 


clg 

semilogx(tpt,ss,‘rX‘) 

xlabel(’log time (sec)’),ylabel(’load (gm)'),grid 
title(’Proof Tested Failure Times’ ) 

pause 


clg 

yl=yl./n; 

semilogx(xtl,yl, y+’) 

xlabel(’log time (sec)’),grid 
title(’Post-Proof Test Safety Zone Histogram’ ) 


% ***e** SAVE DATA ttt 


% Histogram Data 
save XTl.m xt1 -ascii 
save Yl.m y1 -ascii 


* Proof Test Data 
save SS.m ss -ascii 
save TPT.m tpt -ascii 


%* Non-Proof Test Data 
save S.ms -ascii 

save TNPT.m tnpt -ascii 
format short e 


B. PARAMETERS IN POST-PROOF TEST SIMULATION 

A total of 2000 fiber strengths were simulated in each of 
50 runs. A total of nine cases were looked at. For each case, 
the following parameters were used: 
a=5, B=20, L, = .8, bo = .8, Lz = -.8, p = 40, teehat = 1 
t> = 50000, dwell time = 1. 
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7, proof load (dS) = 3 
7, proof load (dS) = 5 
7, proof load (dS) = 7 


10, 
10, 
10, 
13, 
13, 
13, 


proof 
proof 
proof 
proof 
proof 
proof 
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